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1/INTRODUCTION 


This monograph discusses elastic longitudinal models of space launch vehicles and 
methods used to determine the dynamic characteristics needed for analysis of vehicle 
response and stability. The types of vehicles considered are multistage liquid and 
solid propellant vehicles and vehicles with clustered tanks. Longitudinal dynamics have 
been uncoupled from lateral and torsional dynamics in the frequency range of interest 
for most launch vehicle systems studied to date; however, such coupling can be impor- 
tant for spacecraft as well as for launch vehicles with the proper characteristics. 

The basic method of analysis of longitudinal dynamics problems is to formulate 
a mathematical model for the vehicle, establishing mass and stiffness model properties 
from corresponding vehicle properties. Normal modes of the model are then deter- 
mined, and the response to a particular excitation is found from the normal mode and 
excitation properties. This method closely parallels that for lateral dynamics. 

Since analyses are made using properties of the model, the results are only as 
good as the model. The model requirements depend on the particular problem being 
studied and the desired results. Simplifying assumptions are usually made in ideali- 
zing the vehicle, so that only significant dynamic phenomena are represented; usually 
only the lowest frequency modes of the vehicle are of interest. To a great extent, 
judgment requiring some experience and understanding of the vehicle structural com- 
ponents must be used in developing a simple but adequate vehicle model. 

The most significant masses involved in launch vehicle longitudinal dynamics are 
the propellant masses; other important masses related to the vehicle lower frequency 
modes are usually the remaining heavier masses such as payload and engines. The 
significant structural elements are propellant tanks, interstage adapters, and tank 
bulkheads. Vehicles are usually symmetric about their longitudinal axes, and this 
permits the use of what is effectively a one-dimensional model. 

The propellant mass varies slowly with time in the particular stage that is burn- 
ing. This permits the assumption that the mass is constant at any particular instant 
of time. Analysis on this basis is usually referred to as the time-slice approach and 
is adequate for most problems . 

The most significant stiffnesses are usually related to the propellant tanks, in- 
terstage structures, and engine supports. Stiffnesses related to propellant tanks vary 
with time since the basic configuration of the vehicle is changing as propellants are 
consumed. 

For liquid propellant launch vehicles with thin-skinned cylindrical tanks, lumped- 
mass models are adequate for determination of the lowest vehicle modal proper- 
ties. For vehicles with more complex tank geometries, more general methods such 
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as the use of finite shell elements to construct a vehicle model may be required. 
Models of solid -propellant rockets may be obtained by either continuous representation 
of the core and casing or by making a simple model that takes into account the mass 
and shear stiffness of the solid propellant core. 

The problems to be analyzed using longitudinal models include response to engine 
ignition and shutdown transients, release of a vehicle from its launcher, and the 
stability of engine-structure coupled oscillations. A second stability problem 
peculiar to liquid-propellant vehicles with thin-skinned tanks involves coupling between 
structural motions and the tank pressure regulation system. Each of these problems 
may require different emphasis on details of the model used. 
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2 /state of the art 


The analysis of liquid- and solid-propellant vehicle longitudinal dynamics may be dis- 
cussed separately since the modeling problems are of a different nature for the two 
types of vehicles. 

2. 1 LIQUID-PROPELLANT VEHICLES 

The unique problem in liquid-propellant vehicles is representation of liquid-propeUant 
tanks. The longitudinal dynamics of elastic tanks containing liquid propellants have 
been analyzed in considerable detail; for example, Reference 1. The types of tanks 
that have generally been considered are axisymmetric . Effects that have been con- 
sidered include tank initial stresses, liquid surface waves (sloshing), tank surface 
motions that are not symmetric, compressibility of ullage gas, elastic effects of 
bulkheads, and orthotropic tank construction. Analyses including many of these ef- 
fects are very complex; as a result, use of such analyses in developing a vehicle 
model has been limited. 

The tank models used extensively to study vehicle longitudinal dynamics have 
generally included a simplified representation of the tank longitudinal dynamics. Re- 
ference 2, for example, presents various single -degree-of -freedom models for 
longitudinal vibration analysis of liquid-propellant cylindrical tanks. These single- 
degree-of -freedom models were derived based on the assumption that the liquid has 
a longitudinal acceleration (relative to the tank bottom) that is constant from the liquid 
surface to the bottom of the tank. Liquid motion is due to tank wall and tank bottom 
radial and longitudinal deformations. These deformations per unit inertial force of 
the liquid provide information that can be used to develop a model for the tank. Re- 
ference 3 uses a similar tank model to study the coupling between vehicle longitudinal 
oscillations and propellant tank pressure regulation at launch for an Atlas vehicle. 
Reference 4 indicates how to include the effective stiffnesses of bulkheads forming 
tank bottoms for cylindrical tanks. Extension of the single-degree-of-freedom spring- 
mass model to include higher degrees of freedom is also discussed in Reference 3. 

Simple models for other commonly used tank shapes , such as ellipsoidal tanks 
and tanks with conical sections, may be derived in a manner similar to models derived 
for cylindrical tanks, but the problem is somewhat more difficult. 

A recently developed method of liquid tank representation is presented in Refer- 
ence 5. A finite element representation of tank and other vehicle components is used 
to develop the vehicle stiffness matrix; the vehicle is divided into axisymmetric shell 
components. Fluid motions are assumed to be consistent with the shell component 
distortions while fluid sloshing effects are neglected. The vehicle stiffness 
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and mass matrices are obtained by superposition of matrices of individual elements. 
Only axisymmetric shell motions are considered. Each tank may be composed of el- 
lipsoidal, conical, and cylindrical elements. This model may include higher modes 
of the tank, bulkhead effects, and general shapes of tanks. 

Effects not included in the simpler spring-mass models or the model of Reference 
5 are tank motions that are not axisymmetric, i.e. , motions that vary around the tank 
circumference and the effects of fluid surface motion. 

The vehicle mass and stiffness matrices that result from use of the tank models 
discussed may be operated on using standard techniques to determine the vehicle 
natural frequency and natural mode properties to be used in longitudinal response 
analyses. 

The stiffness of a local structure may be evaluated by use of finite elements to 
develop stiffness matrices by methods such as that presented in Reference 6. Local 
structure stiffness should also be verified by tests where possible. 

2.2 SOLID-PROPELLANT VEHICLES 

Solid-propellant rockets have generally been represented by continuous models J Re- 
ference 7 provides a detailed analysis of the structural dynamics of a solid rocket. 

It includes a discussion of propellant stress-strain properties and analysis of several 
problems using a continuous representation. The effective longitudinal stiffness of 
the solid rocket depends primarily on the solid-propellant longitudinal shear stiffness 
since inertial forces acting on the core material are transmitted to the rocket casing 
through shear stresses developed in the core material. The rocket casing may be 
relatively stiff longitudinally in comparison to the core; in this instance the simplify- 
ing assumption of a rigid casing might be made. 

The remainder of the vehicle, adapters, payload, etc. , may be modeled in the 
same manner as discussed for the liquid-propellant vehicle. 
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3/MODEL requirements and recommended procedures 


The problem of determining the reeponse of a structure to a particular forcing hinc- 
tion consists of formulating the equations of motion for the structure and obtaining a 
solution to these equations. The solution is usually easier to obtain if the equations 
of motion are written in terms of the vehicle normal mode coordinates. The use of 
matrix notation also permits simplified discussion of the general analytical problem 
and provides a compact format for organizing data for computations. 

Natural modes and frequencies may be determined from equations of motion of 
a system with external forces equal to zero. The equations of motion may be derived ^ 
from considerations of dynamic equilibrium or by use of energy concepts and LaGrange s 
equation. As a simple example of the use of equilibrium to derive the equations of 
motion, consider the two-degree-of-freedom system of Figure 1. The symbols m and 
K represent masses and spring rates, x represents displacement, and subscripts 1, 

2, and 3 refer to nodes of the structure. Masses are lumped at these nodes, and 
massless springs connect the nodes in this particular example. 



Figure 1. Two-Degree-of-Freedom System 

From a consideration of equilibrium of inertial forces and spring forces un 
free vibration, 


-(Xi-x 2 )K 12 -m 1 x x = 0 

(*! * * 2 > K 12 ' K 23 X 2 ' m 2 H = 0 
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or, in matrix notation, after rearranging terms. 


m J 0 


m_ 


2 J 


K 


12 


-K 


-K 


12 


12 


K 12 + K 23 


= 0 


or simply, 

[M] {x } + [K] { x } = 0 


(3) 


(4) 


where the matrix containing the mass terms is the mass matrix and the matrix con- 
taining the spring or stiffness terms is the stiffness matrix. 


LaG range's equation for determining the equation of motion is 


c>_ d(KE) 
<H dq 


+ *<m 

dq n 



0 


(5) 


where KE is the kinetic energy, PE is the potential or strain energy, W is the work 
done by external forces acting on the system, t is time, and q n is a generalized co- 
ordinate of the system. For the free vibration of the system in Figure 1, 


Ql = x x 

(6) 

q 2 = X 2 

(7) 

KE = }( mi i* ♦ m 2 */) 

(8) 

PE 

(9) 

W = 0 

(10) 


Substitution of the expressions for KE, PE, and W into Equation 5 and letting q = x. 
yields one equation of motion; letting q n = X2 yields the second equation of motion. 
These are the same equations of motion derived from the consideration of equilibrium. 

With the assumption of simple harmonic motion, i.e.. 



where x Q is an arbitrary amplitude and to is a natural frequency, Equation 4 becomes 


-CO 2 [M] {x} + [K] [x] = 0 


( 12 ) 


or 


to = [KT 1 Cm] M 

CO 

= [D]{x] (13) 

With the matrix equation in the above form, iterative methods (see Reference 8) 
may be used to determine the lowest and successively higher natural modes and 
corresponding natural frequencies. For a particular natural frequency co n , a set of x 
displacements which satisfies Equation 13 is the natural mode, [0 n L corresponding 
to cu n , i.e. , at co = co n , [x] = {(^} . [D] is called the dynamic matrix. 

The response of a vehicle to a particular forcing function may be found by sum- 
ming the responses of its natural modes. The equations of motion including external 
forces become 


[M] fx] + [K]{x} = {F} 


(14) 


The displacements {x} may be written as the summation of the responses in the indi- 
vidual modes: 


n=m 


x i = E *ni% 


n=l 

or, in matrix notation, 

{x} = [0] {q} 

where 


[«] = [*! : «ji- 


(15) 


(16) 


f \ 


{q} 




(17) 


Substituting Equation 16 into 14, 

[M] [0] [q] + [Kl [0 j fq] = { F] 


(18) 
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rp 

Premultiplying each term in Equation 18 by [0] 1 yields 

[<S] T [m][ 0 ] {q} +C0] T [K][0] {q} - [0 ] T {F} (19) 

Because of the orthogonality of natural modes with respect to the mass and stiff- 
ness matrices, i.e., 

{0 n l T [M]{<g = fo n } T [K]{0 s } 

= 0 for n / s (20) 

the equations of motion in terms of the natural modes yield one uncoupled equation for 
each mode of the form 

[0 n ] T [M] {0 n } + [0 n } T [K] {0 n } q n = {0 n } T {F} (21) 

rp 

The terms {0 n }^ [M] {0 n } and {0 n 3 [K] {0 n } are usually referred to as the generalized 

mass and stiffness of the n^ 1 node and expressed as 

\ = {0 n ] T [M] [0 n ] (22) 

= ^nJ T [K] ^n^ (23) 

Noting that 


2 

oo 

n 



(24) 


is the square of the natural frequency of the rfi 1 mode and letting the generalized force 
of the rft 1 mode be defined as 


Q n (t) = {0 n } T [F} 


Equation 21 becomes 


% 


+ co 


% 




(25) 


(26) 


The modal amplitude (generalized coordinate) q n may be found for each mode as 
a function of time from Equation 26, and then the total response may be found from 
Equation 16. 
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Then, to determine the response of a structure using modal properties, a model 
that yields mass and stiffness matrices for determining modal properties is desired. 
The selection of the specific model must be made in a manner such that modal pro- 
perties are accurately defined where forces are applied and such that resulting modes 
give accurate values of the generalized mass. The model must be defined in enough 
detail to provide accurate response information in the area of interest on the vehicle. 

3. 1 LIQUID-PROPELLANT VEHICLE 

A vehicle may be made up of several stages, with the major components by weight 
being propellant tanks, engines, payload, and fairings. Propellant tanks represent 
important portions of the liquid-propellant vehicle model since propellant weight is a 
large part of vehicle weight throughout most of the flight. Tank models are parti- 
cularly important at liftoff when the tanks are full and for vehicles with only a few 
stages. 

Tank models may be of the lumped parameter (spring-mass) type, which are 
often adequate to represent the dynamic characteristics of certain types of tanks, 
particularly cylinders. Continuous representations of the tanks may also be used. 

The method of Reference 5 uses sets of assumed displacement functions for tank 
elements in a manner that takes into account tank shapes made up of cylindrical, 
ellipsiodal, and conical elements. 

The models for structural elements not containing liquid may be represented by 
simply the effective longitudinal stiffnesses of the structural element. This stiffness 
may be determined by tests or by detailed analysis using methods for developing 
multidimensional stiffness matrices of structures, such as the methods discussed in 
Reference 6. 

In the following paragraphs, spring-mass models for various forms of cylindrical 
tanks and models based on References 2 and 3 are discussed. 

3.1.1 SELECTION OF MASSES. The number of degrees of freedom of a model 
determines the theoretical number of modes and frequencies. However, for a model 
with a large number of degrees of freedom, only a limited number of modes and fre- 
quencies may be determined accurately, and only the characteristics of the lowest 
modes are generally of interest. As mentioned above, the major portion of vehicle 
weight is propellant weight throughout much of the vehicle flight. The model of each 
propellant tank must represent the significant longitudinal modes of the tank in the 
frequency range of interest for the vehicle. This frequency range is generally from 
5 to 30 Hz for liquid-propellant launch vehicles. When lower stage liquid tanks are 
full, several longitudinal degrees of freedom of the propellant mass might need to be 
included in a model to represent the lowest four or five elastic modes of the vehicle. 
As propellant tanks become more nearly empty, a single -longitudinal -degree-of- 
freedom representation of a propellant mass may be all that is required to compute 
the lowest four or five vehicle modes. 


9 



The details of the representation of other portions of the vehicle are not usually 
as critical as the propellant representation. Uncoupled longitudinal frequencies of 
engines, payloads, and fairings maybe well above the lower vehicle frequencies, and 
therefore these components might be represented by a single lumped mass, or in 
more detail if desired. In a vehicle with a large number of stages the interstage 
stiffnesses may be relatively flexible and have a more dominant effect than upper 
stage tank stiffnesses on the vehicle lower modes and frequencies. In such instances 
upper stages might be represented by lumped masses connected by the interstage 
stiffnesses. In other words portions of the vehicle which are relatively rigid (have 
high uncoupled natural frequencies) might be considered as lumped masses, whereas 
portions of the vehicle which have low uncoupled natural frequencies and large masses 
must be modeled in enough detail to represent these frequencies. 

3. 1.2 LIQUID PROPELLANT AND TANK REPRESENTATION: SPRING-MASS 
MODELS. In the discussion below, simple single-degree-of-freedom spring-mass 
model representations are developed for cylindrical tanks. The flexibility of an el- 
liptical tank bottom is included based on the results of Reference 4. Effects of 
stringers and buckled skin are considered as variations of the basic model. A model 
for small ullage is developed for use in analysis of tank pressure regulation systems. 
Models for several degrees of freedom and for cone-cylinder tanks are also con- 
sidered. These models are based on models developed in References 2 and 3. 

The analysis of a coupled elastic tank and propellant mass, even for a highly 
simplified case, becomes a very complicated eigenvalue problem. Specifically, a 
solution must be obtained that satisfies the differential equations for the liquid and 
the elastic shell, as well as appropriate boundary conditions at the tank walls, the 
tank bottom, and the liquid free surface. A rigorous analysis of this type is reported 
in Reference 1. The results yield the natural frequencies for the tank and propellant. 

A number of other analyses have also been attempted for the coupled liquid and 
elastic container. Most of these generate a great deal of mathematical analysis that 
is of very little use in defining an analytical model for the tank. It is apparent, then, 
that other, more simplified techniques must be used. The continuous analysis can then 
be used as a check on Hie dynamic characteristics of the simplified representation. 

A major building block in the development of a longitudinal structural model is a 
lumped parameter model for each propellant tank. In general, such tank models have 
been restricted to a single (predominant) mode of the coupled elastic shell and propel- 
lant mass. However, there is reason to believe that a single-mode model is not ade- 
quate in all cases. Furthermore, the development of a multimode model for the pro- 
pellant tank is feasible. 

As a basis for developing a tank model, first it is necessary to establish a set of 
approximate equations for the shell and the liquid. For example, consider the tank 
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shown in Figure 2 # filled to a height L with a 
nonviscous, incompressible liquid, where the 
axial acceleration of the fluid is x. If it is as- 
sumed that the tank shell is thin, that any ef- 
fects of preloading can be neglected, and that 
the loading is axially symmetric, then the fol- 
lowing equations can be developed (Reference 
10 ). 

2 4 2 r 

12(1-10 dX a 


+ — N = p 
a X F 


r = a 


% 

dX 


- v v 

— ; — N„„ W 

Eh X a r 


(27) 



Figure 2. Cylindrical Tank and 
Liquid 


(28) 


where 


w^ = the radial displacement of a shell element 

the longitudinal displacement of a shell element 
the axial force/unit circumferential length 
the liquid pressure 
Young T s modulus 
Poisson’s ratio 
a = the shell radius 


u = 
X 

N x = 

p = 
E = 
v = 


h = the shell thickness 

Equations 27 and 28 are essentially the same as those defined for a cylindrical 
shell (Reference 9) and a static loading condition. If, in addition, the shell is very 
d^w 

thin, the — term contributes very little to the gross deformation of the shell. 

dX 4 

Neglecting this term, Equation 27 becomes 
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(29) 


w 

r 



ay 


Eh N X 


For a nonviscous, incompressible liquid, the fluid velocities are defined in terms 
of a velocity potential <t> by Laplace's equation, i.e. , 


d 2 g> 1 d<t> + j_ d 2 $> 

„ 2 r dr 2 ,2 2 

dr r 30 33T 

where 

d$ jdr 
dr dt 

1 3<1> _ d0 
~ r d 0 ~ F dt 

d<t> _ dx 

" dX dT* 


(30) 


When all forces applied to the fluid are axially symmetric, <t> is not a function of the 
coordinate 0 and Equation 30 reduces to 



1 d<E 
r dr 



(31) 


For small fluid velocities, the total fluid pressure is given by 


P 



x (X- L) 


(32) 


Equations 28, 29, 31, and 32 are used along with the appropriate boundary con- 
ditions of the liquid surface, tank wall, and tank bottom to form the basis for the 
analysis given in Reference 10. In this analysis, the effects of shell inertia for a 
thin shell are neglected in comparison to the liquid inertia, such that Equations 28 and 
29 also apply to the dynamic condition. 


In order to further simplify the development of a lumped parameter model, the 
velocity potential term in Equation 32 is usually neglected; i.e. , the following approxi- 
mation is made: 


p = p x (L - X) 


(33) 
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Equations 28, 29, and 33 are then used to develop the tank model. This neglects 
all contributions due to the liquid velocities, as defined by the velocity potential 4>, and 
the liquid-free surface waves. 

Equations 28 , 29 , and 33 can also be used as a basis for developing a multimode 
tank model by dividing the tank into two or more sections. This approach can account 
for the higher tank modes as well as variations in tank geometry and skin thickness. 

The propellant tank most commonly encountered is usually an internally pres- 
surized circular cylinder with a rounded bottom. In many cases, the bottom is ap- 
proximately elliptical in shape and can be either concave upward or downward. 


In this section, a basic single-mode tank model is developed for the case of a 
large tank ullage volume where the pressure essentially remains constant. This model 
is then modified to account for the individual effects of stringer reinforcement and a 


buckled skin with stringer reinforcement. 

3. 1.2. 1 Basic Single-Mass Model. The tank shown 
in Figure 3 is partially filled with an incompressible 
liquid to a depth L and the tank ullage is pressurized 
with a gas to some constant pressure. The tank is an 
elastic circular cylinder with a radius a and a con- 
stant thickness h. The tank bottom, also elastic, has 
an elliptical shape, where the major semi-axis is a 
and the minor semi-axis is b. The bottom can be 
either concave upward or downward. In addition, it 
is assumed that all loading on the tank and liquid is 
axisymmetric. 



.. Figure 3. Partially Filled 

For an axial acceleration of x, the liquid pres- Tank 

sure acting upon the wetted portion of the tank is 

given by Equation 33. Hoop stresses are produced in the tank skin due to this pres- 
sure. These stresses result in radial and longitudinal displacements for each ele- 
ment of the tank with respect to the bottom. The radial displacement, which is 
proportional to the pressure, is shown in Figure 4a. 


When the tank is assembled with a vehicle, an acceleration of the masses above 
the tank results in a longitudinal force acting on the top of the tank. This force pro- 
duces an additional axial stress in the tank skin and displaces the tank and the liquid 
center of gravity. 


A flexibility matrix can be developed relating displacements at X = L and center 
of gravity of the tank to a force at X = L and an effective ’’liquid inertial force” acting 
at the center of gravity of the liquid when the bottom of the tank is fixed as shown in 
Figure 4. A stiffness matrix and a corresponding model for the tank can then be 
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Figure 4. Tank Strains 


developed. This approach yields a single-mass tank model similar to that obtained in 
References 10 and 2. 


The equivalent single-mass model is indicated in 
F igure 5 . Nodes 1 and 2 of the lumped parameter model 
are used to signify* respectively, the tank at X = L 
and the liquid center of gravity. The displacements 
of these points relative to the tank bottom (i.e. , 
relative to X = 0 in Figure 3) are denoted by x^ and 
X£. The force acting at the top of the tank is F^ and 
the effective force acting at the center of gravity is 
F£, the total inertial force of the fluid due to x , i.e. , 



Figure 5. Equivalent Tank 
Model 


F = Vpx = m x 

Ck 


(34) 


where V, p, x, and m are the tank volume, fluid mass density, fluid acceleration, and 
the fluid mass, respectively. The liquid volume V for the cylindrical tank with an el- 
liptical bottom is given by 


V 


n a LD 


(35) 


where 


D 


1.4^. 

3 L 


(36) 


Therefore 

F = fra 2 LDpx 


(37) 
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( 38 ) 


^ 2 

n a LD 

It should be noted that when the bulkhead In Figure 3 is inverted (the dashed line) , b 
takes on a negative value but the above expressions remain valid. 

The relationship between displacements and the liquid inertial force are developed 
first. The longitudinal and circumferential stresses and strains in the thin tank skin 
due to the fluid pressure, p, are (to a close approximation) given by 


} (39) 


where a x and oy are the tank longitudinal and hoop stresses, e x and Cy are the longi- 
tudinal and hoop strains, v is Poisson's ratio, and E is Young's modulus. Positive 
stresses are defined to be compressive in the longitudinal direction and tensile in the 
hoop direction; positive displacements of the tank top and liquid center of gravity are 
downward. 

The longitudinal displacement of the tank at the liquid surface due to the fluid pres- 
sure, p, is 

L 

*12 * / C x dX («) 

o 

where e x is given by Equation 39. Thus, substituting Equation 33 into 39, 

L •• * 2 •• 

r ivapx(L-X)dX iv a L px 

X 12 = J Eh = — 2Eh < 41 > 

o 

Substituting in the expression for p x from Equation 38 yields 
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X 12 " 2tr DE h a 


(42) 
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The tank wall radial displacement, w , due to p (Equation 33) is given by 

r 

2 2 

pa pa (L -X) x 

w = a f — — = — — =-r - 

r y Eh Eh 


( 43 ) 


The displacement of the center of gravity due to w r and the tank bottom stiffness is 
then approximated by 
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22 


2 L 

3 V 



277 aw 
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dX + 



o - T 2 F „ 

2 pxaL t 2 

3 DEh + K~ 
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(44) 


where Kg is the tank lower bulkhead stiffness (the appropriate values for Kg are de- 
fined parametrically later). Substituting Equation 38 for p x into Equation 34 yields 


x 


22 


2LF„ F . 

2 2 _ 

3rraD 2 Eh K b 


(45) 


Considering an axial load acting on the tank, as shown in Figure 4b, F^ is an axial 
load in the tank skin resulting in a load per unit circumferential length of the tank skin 
given by 



The stresses and strains in the tank due to F^ are 


(46) 


a = - — r 
x 2?rah 

a - 0 
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x 2iraEh 
v F 

1 

y 2fraEh 

while the displacements of the tank at X = L and center of gravity due to F^ are 


(47) 
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L 

= / 


LF. 


dX = 


(48) 


27raE h 


16 



X 


21 



2 JV 1 
2ira € dX • — 

y V 


vLF i 

27raDEh 


The axial stiffness of the tank skin is 


( 49 ) 


_ 2traE h (50) 

K = L 

Then the relationships between forces and displacements from Equations 42, 45, 48, 
and 49 are 
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1 _ 
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= = Ji_ 

12 ~ F DK 

= 111 = Ji_ 

d 21 F " DK 


d 


22 



4 1 

3 D 2 K K B 


( 51 ) 


or, in matrix notation, 




[d] | F | 


( 52 ) 


The [d] matrix in Equation 51 can be inverted to obtain a stiffness matrix for the tank. 
The springs of the equivalent model for the tank can be obtained from this matrix. 
From Equation 51, Fj and F 2 are 
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-3l/ DK 


( 54 ) 


- _ - 3 1/ DK 3 D 2 K 

2 d 1 2 o r» 2 K + yl 2 2K X 2 

4-3l'+3D — 4-31/ + 3D 


K 


B 


K 


B 


or, in matrix notation, 
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K 21 K 22 


Therefore, 
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(56) 


The springs of the equivalent tank model shown in Figure 5, determined from the 
coefficients of the stiffness matrix (as indicated in Reference 6), are 
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Thus, 
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(Contd) 


K_ = 


3D(D- l/)K 

. 2 2 K 

4-31/ + 3D — 

k b 


K = 


4- 31/ D + 3D 2 7T- )K 

i 

4 - 3 l/ 2 + 3 D 2 — - 

k b 


(59) 


3. 1.2.2 Spring Rate for Elliptical Tank Bottom. In Reference 4, spring constants 
are presented for ellipsoidal tank bottoms , as determined from an analysis using 
linear membrane theory. Graphs from which the spring constants can be obtained 
are reproduced here. 

Consider, once again, the tank shown in Figure 3. From Reference 4 die tank 
bottom spring rate is given by: 


K = Eh 2T<3e ' 2t)2 T < 60 > 

9 1 H (f, v) - 2 eG (f, l/) + e F (f, v ) J 

where 

E = Young's modulus 

h = bulkhead thickness 

B 

a = radius of cylinder 

b = bulkhead semiminor axis 

L = height of liquid in cylindrical portion of tank 

e = L/a 

f = b/a 

v = Poisson's ratio 

The functions H(f,i/), G(,P), and F(f,l/) are given in Figures 6, 7, and 8. These 
functions are derived in Reference 4 . 

When the tank bottom is inverted as indicated by the dashed lines in Figure 3 the 
expression for the spring rate becomes 
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F<*. v) 


( 61 ) 


K' = E h 
B B 


2 7T (3 e + 2 f) 


9 |h (f, v) + 2 eG (f, u) + e 2 F (f, i/)J 



Figure 6. Variation of F(f,i>») with Depth- 
to-Radius Ratio for Various 
Values of Poisson f s Ratio 


Figure 7. Variation of G(f,v) with Depth - 
to-Radius Ratio for Various 
Values of Poisson’s Ratio 


3. 1.2.3 Tank with Stringers and Buckled 
Skin. A model cylinder with skin and 
stringer construction is shown in Figure 
9. If the stringers have negligible radial 
stiffness and the skin is unbuckled and 
can displace longitudinally and radially, 
the tank models described in the preceding 
subsections can be combined in parallel 
with a spring K g , where 

_ (AE) . 

K = stringers (62) 

s L 

is the total axial stiffness of the stringers 
in the cross-section. 

When the skin is partially buckled in 
the axial direction, the axial stiffness of 
the tank skin is 



Figure 8. Variation of H(f,i/) with Depth- 
to-Radius Ratio for Various 
Values of Poisson’s Ratio 
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Figure 9. Skin-Stringer Cylinder Tank Model 


and the flexibility coefficient d-Q in Equation 51 becomes 


11 


K 


Bu 


+ K 


(64) 


The center-of-gravity displacement due to the axial load Fg u carried by the buckled 
skin is then (neglecting the bulkhead and liquid compressibility effects) 


x 


21 


v F 


Bu 


K 


(65) 


where K is the axial spring constant for the skin when unbuckled. Equation 50. That 
portion of the load carried by the buckled skin is given by 


K F 
Bu 1 


Bu 


K 


Bu 


+ K 


( 66 ) 


where F^ is the total axial load carried by the skin and stringers. Then 
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and 


VK F 
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K (K + K ) 
Bu s' 


(67) 
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( 68 ) 


The center-of-gravity displacement due to acceleration of the liquid in the tank is still 
(neglecting bulkhead and compressibility effects) 


d 


22 


4 

3K 


(69) 
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Consequently, the flexibility matrix becomes 


[d] = 


K + K 
Bu s 


uK 


Bu 


K (K + K ) 
' Bu s' 
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K (K + K ) 
Bu s 
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( 70 ) 


After inverting the [d] matrix and solving for the spring rates, the values for the 
model in Figure 10 are 
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(71) 



Figure 10. Skin- Stringer 
Model with 
Buckled Skin 


3. 1.2.4 Effects of Small Ullage Volumes. When the tank ullage volume is small, 
any longitudinal oscillation of the structure can produce a corresponding oscillation of 
significant magnitude in the ullage pressure. This can be the result of a longitudinal 
force transient such as thrust buildup or launcher release. 

The change in ullage pressure is due to a change in ullage volume and ullage gas 
weight (if a pressure regulator is involved). In this section, a single-mode tank model 
is developed for such a case. The effect of gas compressibility is included in the model 
spring rates while the pressure change associated with a change in gas weight is used 
as the forcing function for the tank. 

In order to formulate an analytical representation for the perturbation changes in 
the tank ullage, the following assumptions are made. 


22 



a. The gas, or mixture of gases, in the tank ullage behaves in a quasi- static manner, 

b. The gas process is adiabatic. 


c. The perturbation variables are small. 

The equation of state for a gas (or a mixture of gases) at low pressure is given to a 
close approximation by the equation below: 

PV = WR T ( 72 ) 

m 

where P is pressure, V is volume, W is weight, T is temperature, and R m is the gas 
constant for the mixture. In a strict sense, Equation 72 applies to a condition of ther- 
modynamic equilibrium. It also applies (approximately) for a perturbation condition 
when the changes in these variables are slow enough, and small enough, such that the 
gas is at all times close to thermodynamic equilibrium. The gas is then said to be- 
have in a quasi-static manner. 

For an additional relationship between the variables in the tank ullage, it is con- 
venient to assume that the process is adiabatic. In such case, 

(Vi) 

T oc p y m/ (73) 

where y m is the adiabatic exponent for the gas. In general, a polytropic exponent 
would be used for Equation 73. However, in the frequency range of interest it can 
be assumed that there is little chance for heat transfer to the gas so that the poly- 
tropic exponent is close to the adiabatic value. 

Equations 72 and 73 can now be combined to form the following equation: 


P = 


const. 



(74) 


A given tank may contain a mixture of gases. The total weight is then given by 


W = W + W + ... 

1 tt 


(75) 
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and the effective gas constant, R m , is defined by 


R 

m 


W R + W R + ... 
1 1 2 2 

W 


( 76 ) 


where Rp R 2 , etc. , are the gas constants for the individual gases. 

For small perturbations during flight, the change in tank pressure is defined by 
the first terms of a Taylor's series, where the partial derivatives are evaluated at 
the steady- state (or equilibrium) condition, i.e., 


AP 


d P 
SW 


AW 


BP 

av 


AV 


J s. s. 


s. s. 


(77) 


where AW is the change in ullage gas weight and AV the change in ullage volume. 
For the steady- state condition, 


P = 



(78) 


The partial derivatives in Equation 77 can be evaluated by using Equation 78 to define 
the proportionality constant, so that 


AP = p = K — - K v 
a w s v 

where s is the Laplace operator, and 


(79) 


K 
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K 
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yP 
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(80) 


• • 

where Rg is the gas constant for the pressurizing gas and where w = AW, v = AV, and 
the bars denote steady-state quantities. 
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Figure 11 shows a thin skin, cylindrical tank shell with an elliptical bulkhead for 
the bottom. The top of the tank is assumed to be rigid and the tank is nearly filled 
with liquid. Forces Fj acting on the top of the tank, an effective inertial force, Fg, 
acting at the liquid center of gravity, and an ullage pressure change, p^, acting in the 
ullage space are also shown. The pressure change p^ represents only that portion of 
the total pressure change in the ullage that can be attributed to a change in gas weight. 
Hie compressibility effect, due to a change in ullage volume, is included in the model 
spring rates. Therefore, from Equation 79, is defined by 


or 


sp = K w 
a w 


p a = 


p + K v 
a v 


(81) 


(82) 



Figure 11. Tank with Small Ullage Volume 
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If the displacements of the tank top and liquid center of gravity are and x 2 , 
and the change in ullage volume is v a , then the relationship between forces, pressure, 
displacement, and change in ullage volume is given by the matrix equations 
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: — is the compliance 
K v 

1 

d 

d 

11 

12 

lv 

1 

d 

d 

21 

22 

2v 

1 

d 

d 

vl 

v2 

w 


= Cd] 


26 


1 




1 

1 


K 


DK 


1 

1 

K 





1 

1 

1 

2 ” “ 2 (' -l) 





1 

1 

DK 

V 


4 + 

l 

1 

1 

2 

n a 

DK 


| 3 D 2 K 


1 

1 

+ k b 

— 

— 

| 

1 2 it a 2 /l 

| ' 

-7) 

1 

1 

1 

2 4/5 „ \ 

2 " a \2~ 2V ) 



1 DK 


1 

K 

2 it a 2 (v 

-t) 

1 

i 


1 


K 


1 

1 it a 2 


1 

1 

1 

2 4 
. it a 



! 




( 85 ) 


The compliance, C, could be defined to include flexibility of the tank upper bulkhead. 
As the ullage volume increases, C becomes infinite and the model reduces to the basic 

tank model. 

The [d] matrix can be inverted to obtain a stiffness matrix: 
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A tank model where the ullage gas compressibility is included in the spring rates 
where the applied forces are due to p^ is shown in Figure 12* 
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Figure 12. Tank Model for Small 
Ullage Volume 


The tank spring rates, defined in terms 
of the stiffness matrix, are 
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3. 1.2. 5 Multimass Models. In most analyses of current vehicles, the structural 
representations used have been based upon the tank model shown in Figure 5. This 
model provides an approximation for the first (predominant) mode of the coupled pro- 
pellant and elastic tank. However, it is known that higher modes exist in the tank and 
a limited number of these may have a significant energy content, compared to the first, 
and they do provide a resonant condition at discrete higher frequencies. This is 
shown in Reference 10. 


This implies one of two conditions. Either the higher tank modes are unimportant 
and their omission does not affect the overall structural modes, or certain structural 
modes are inaccurate because of this omission. The actual condition, of course, de- 
pends upon many factors such as the type of structure, the propellant level, etc. 

There has been very little opportunity to determine the accuracy of predicted ve- 
hicle modes using current models. It is believed that, in most cases, the first mode 
predictions are adequate and that any errors are incurred in the higher modes. This 
should be particularly true when the propellant mass is large and represents a major 
part of the vehicle mass. 


It is desirable to further explore the importance or contribution of the higher tank 
modes by developing a multimode tank model and using it in the structural representa- 
tion for the vehicle. The initial results obtained by using this approach are discussed 
in this subsection. 


A discrete model for higher tank modes can be developed in a manner similar to 
that used for the single-mass representation. Consider again the cylindrical tank with 
an elliptical bulkhead as shown in Figure 13. Assume, for the present, that the bulk- 
head is inelastic (rigid) and that the liquid is incompressible. 
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Let the tank and propellant be divided 
into two elements. A point on the tank at 
the liquid surface, the center of gravity of 
the upper element, and the center of gravity 
of the lower element are denoted by Sta- 
tions 4, 3, and 2, respectively. The pro- 
pellant heights are Lj and L 2 while the tank 
wall thicknesses are h^ and hg. 

The procedure used to develop a two- 
mass model is exactly the same as that used 
for the single-mass model. This procedure 
can be expanded to any number of elements. 


The following quantity is similar to 
that from Equation 36. 


, 2 b 

°1 1 " 3 L 1 


( 88 ) 


while the axial rates at the two tank elements are given by 


2 v a Eh 



2 tt a Eh 



(89) 


Next, consider the individual effects of an axial acceleration, x, on the two fluid 
elements ! Hie resulting radial deflections of the tank are shown in Figure 14a and b 
along with the corresponding inertial forces, F 2 and Fg, actingat the two fluid center- 
of-gravity locations. Finally, an external force, F 4 , is considered at the top of the 
tank. This produces the radial tank deflection shown in Figure 14c. The individual 
axial and radial deflections of the shell can be evaluated using the shell equations. 
This leads to the following influence coefficients . 
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The influence matrix is then given by 
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where for L = L = — , Equations 88 and 89 become 
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For a particular tank, the coefficients of the [d] matrix can be evaluated. The [d] 
matrix can then be inverted to obtain a stiffness matrix. Figure 15 shows a model 
which would correspond to the stiffness matrix given by 
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Figure 15. Two-Mass Tank Model 



This model can be generalized for many degrees of freedom by writing general expres- 
sions for the influence coefficients of a multimass model. Also, the effects of a flexible 
tank bottom and liquid compressibility can be included in the process. 
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3. 1.2.6 Tanks with Conical Sections. A tank may be composed of a combination of 
circular cylinder and conical segments where a change in vehicle diameter is re- 
quired between stages. Figure 16 shows a typical example. 



Figure 16. Vehicle with Cone-Cylinder Tank 

An approximate model for a tank of this type might be determined using a cylin- 
der having the same mass and total axial stiffness. For the tank of Figure 17, the 
axial stiffness is 
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( 96 ) 
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Figure 17. Cone-Cylinder Tank 


cylinder axial stiffness 
2 tt a fc) Eh 
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If m cc is the total mass of the liquid in the tank, and nigg and Kgg are used with the 
model of Figure 5 and Equations 56 and 57 , 
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K 

3 4-3V, 2 C ° 

where the model is as indicated in Figure 18, 

This model yields only approximately the correct 
frequency; use of such a model in a complete vehicle 
indicates that its frequency is too high. For a 10- 
degree cone with L^/L^ of about 0.4 for a first-stage 
LC >2 tank, the resulting vehicle first-mode frequency 
is about 20 percent too high in the first mode when 
compared with flight data. As the liquid surface 

drops in the cone and approaches the top of the cylinder , the model frequency appears 
to be more accurate. The apparent error of this model seems to be related to the 
geometry of the cone rather than its stiffness. 



Figure 18. 


Equivalent 
Cylinder Model 


An approximate model for the cone -cylinder tank could also be determined using 
a two-degree -of -freedom model for the tank and assuming that the cone is rigid when 
determining effects of propellant inertial forces. A model of this type is shown in 

Figure 19. 

The same notation is used as in Subsection 3. 1.2.5 except that now Station 3 de- 
notes the center of gravity of the liquid in the conical section of the tank. Now con- 
sider the individual effects of an axial acceleration acting on each of the two fluid 
elements. The resultant radial deformations of the cylindrical portion of the tank are 
shown in Figure 20a and b along with the corresponding inertial forces, F 2 and F 3 , act- 
ing at the two fluid center-of -gravity locations. The effect of an external force, F 4 , is 
indicated in Figure 20c. The axial and radial deflections of the cylinder can be evalu- 
ated using equations for strain developed in a manner similar to that of Subsection 
3. 1.2. 1. This leads to the following approximate influence coefficient matrix: 
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If we want to compute the natural frequency of the tank, the displacement x 4 need 
not be considered: 
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A nontrivial solution of the above equation exists only if the following determinant is 
equal to zero: 
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where I is an identity matrix. 


After substituting expressions for the influence coefficients and making the 
following substitutions, 
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then the determinant becomes 
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Expansion of the determinant leads to the characteristic equation, 


X - b X + c = 0 


(109) 


where 


b = 1 + 3/3 C 


c = TP C 
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Two roots for X result from the characteristic equation. The root corresponding 
to the lowest natural frequency leads to 
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where 


R = 30C Q 

A similar equation may be derived for the "equivalent" cylinder; letting the first 
natural frequency from the equivalent cylinder be 



where 



The two equations given above were evaluated for a ft range of 0 to 0.5 for a tank 
with the following parameters: 

a = 10 ° 

a^ = 60 in. 

= 310 in. 

The results are shown in Figure 21. The two-degree-of -freedom model gives a 
frequency increasingly lower than the equivalent cylinder model as £ increases, i.e. , 
as the liquid level rises higher and higher in the cone. The value of/J = 0.42 corres- 
ponds to about the propellant configuration of an Atlas LO2 tank at liftoff. For this 
value of /3 , the natural frequency given by the equivalent cylinder model is about 17 
percent higher than the value given by the cone-cylinder model. 

The equivalent model for the cone -cylinder tank would be obtained by inverting 
the flexibility matrix of Equation 101 (in a manner similar to that of Subsection 
3. 1.2.5). 



Cx'o = - = FIRST MODE FREQUENCY 

^ pa L L OF CYLINDRICAL SECTION 
b 1 




Figure 21. Natural Frequency of Cone-Cylinder Tank 
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The corresponding spring-mass model (see Figure 19) is 
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3.1.3 LOCAL STRUCTURE EFFECTS. The stiffness of local structure may be 
evaluated in detail using methods for representing complex structure such as that 
presented in Reference 6 or by direct tests after the hardware has been built. Tests 
on local structure are advised where practical since difficulties are often encountered 
in representation of support structure for engines, payloads, and components. 


3.1.4 TEMPERATURE. Extreme temperatures ranging from cryogenic to several 
hundred degrees Fahrenheit may affect the value of the modulus of elasticity of impor- 
tant portions of the vehicle structure. The correction required may affect frequencies 
by several percent at certain times of flight and should be taken into account if this 
degree of accuracy is required. 

3.1.5 EFFECT OF AXIAL LOAD ON STIFFNESS OF SKIN-STRINGER STRUCTURES. 
This subject was covered in Subsection 3. 1.2.3 in the discussion of a tank with strin- 
gers and buckled skin. Any skin and stringer structure has an axial stiffness that is 
flie sum of the stiffnesses of the stringers plus the effective stiffness of partially 
buckled skin between the stringers. 


3.1.6 IMPROVED ANALYTICAL MODEL. The model of Reference 5 uses a finite 
element technique to construct the total launch vehicle stiffness matrix [K] and mass 
matrix [M] by dividing the vehicle structure into axisymmetric shell, fluid, and 
spring-mass elements. 
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The Rayleigh-Ritz technique is used with assumed polynomial displacement 
patterns to determine shell and fluid element mass and stiffness characteristics. The 
stiffness and mass matrices for the complete vehicle are obtained by superposition of 
the stiffness and mass matrices of the individual shell, fluid, and spring-mass 
components . 


In the Rayleigh-Ritz method a stiffness matrix is determined for each element 
from the strain energy of the element. The strain energy of the element can be 
determined from the assumed displacement functions and the stress-strain laws as- 
sumed for the shell element. The strain energy for a shell of revolution due to axi- 
symmetric loading is 


V = 
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where r and r and the radii of curvature of the shell in the meridional and hoop 


directions. N is the initial meridional stress, and 
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where 
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C^4 are orthotropic constants , 

= meridional and hoop load per unit length 
= meridional and hoop moment per unit length 
= meridional and hoop strain 
= meridional and hoop curvature 


Similarly, mass properties for fluid and shell elements may be computed from 
their kinetic energies based on the assumed displacement functions. This method, in 
the current programmed form of Reference 5, has been applied to analysis of a one- 
tenth scale Saturn V model as reported in Reference 11* 

The method appears to be the most general method available for analysis of a 
complete vehicle. It overcomes many of the limitations of the approximations re- 
quired in developing the discrete lumped parameter models discussed earlier. In 
its current form (Reference 5) it will not handle such tank configurations as the cone- 
cylinder tank discussed in Subsection 3. 1.2.6; it does not account for ullage gas 
effects, etc. However, the generality of the basic method would conceivably permit 
the model of Reference 5 to be extended to account for such effects. 

3.2 ADDING COMPONENTS USING MODE SYNTHESIS 

Frequently it is desirable to make a parameter study to determine the effect on ve- 
hicle response resulting from changes in the characteristics of a specific area or 
component, e.g. , a sloshing mass or engine system. Rather than make several 
analyses of the system, changing but a fraction of the parameters each time, the 
vibration characteristics of the system excluding the specific varying parameter may 
be calculated, and then modified by coupling the parameter back in through the mode 
synthesis technique (discussed in References 12 and 13). 

The equations of motion may be written for all components of a system. The 
equations may be written compactly in the uncoupled form for components of a system 
having an arbitrary number of components. Interpreted physically, this can be con- 
sidered to be a set of equations of motion for a group of components which are not 
connected. A set of relationships exists which corresponds physically to connecting 
the components to form a system. These relationships result from requirements 
that displacements at mutual attachment points must be equal. This corresponds to a 


44 



tr ans formation relating the component coordinates to the system coordinates. The 
details of the transformation depend on the system being analyzed. 

This same discussion applies in a similar manner to normal modes of components. 
The motions of one component may be expressed in terms of its normal modes while 
the motions of a second component may be determined in terms of its normal modes. 
The normal modes of the combined system may be determined in terms of the compo- 
nent modes by satisfying the requirements for compatibility of displacements at the 
attachment points between the two components. 

This method would be useful, for instance, for determining the natural modes 
and frequencies of a launch vehicle with several different payload configurations. The 
natural modes of the vehicle could be computed for the vehicle without the payload and 
for the payload alone; the combined system could then be determined in terms of these 
modes for each payload configuration. 

The method offers the advantage of being able to combine analytically determined 
modal data with experimentally determined modal data where available. For large 
structures, modes of individual portions may be determined by test, while tests of the 
entire structure may not be practical. 

3.3 CORRECTING MODEL BASED ON TEST RESULTS 

The final evaluation of analytical techniques is a comparison with experimental data. 
Perfect comparisons are indeed exceptions, since both the analytical model and ex 
perimental model are approximations to some extent. The analytical approximations 
have been discussed. The major experimental approximations are centered around 
suspension system effects and vehicle modifications required to accommodate the 
suspension system. No general rule can be made to obtain better agreement between 
test and analysis. Careful examination of the data and the structure will probably in- 
dicate several areas where the representation is inadequate or does not define the 
test specimen. Possible causes of differences are: 

a. Effects of suspension system on test environment. 

b. Stiffness of joints or trusses. 

c. Assumed planes of symmetry are incorrect. 

d. Effect of large components such as engines. 

e. Experimental modes maybe impure, i.e., not orthogonal. 

f. Effects of moment of inertia. 

g. Nonlinearity. 
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The work of Reference 14 presents a method for obtaining the flexibility matrix 
from experimental mode data. The procedure orthogonalizes the experimental modes, 
using an analytical mass distribution, and then derives the flexibility matrix of the 
structure. This method can be useful if complete and accurate experimental data are 
obtained for a system difficult to model. It can also be used to locate possible dis- 
crepancies between analytical and experimental results. 

3.4 SOLID- PROPELLANT BOOSTERS 

The structural dynamics of solid propellant rockets have generally been treated by 
continuous representations, for instance as in Reference 7. The effective longitudinal 
stiffness of the core depends primarily on the shear stiffness of the solid propellant. 
Inertial forces acting on the core material are transmitted to the rocket casing through 
shear stresses developed in the core material. An element of a solid -propellant rocket 
is shown in Figure 23. 





Figure 23. So lid -Propellant Rocket Element 
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The inertial forces p x acting on the core cause shearing of the core relative to 
the casing and produce an increment of axial force in the casing, over the length 

dX. The shear deformation is indicated under axial acceleration of the element. A 
simple model for the solid-propellant rocket could be obtained by integrating the ef- 
fects of axisymmetric shear deformations of the propellant and axial deformations of 
the casings. The propellant may be relatively more flexible than the casing; the 
casing might then be assumed rigid. 

Reference 7 provides a more detailed analysis of the structural dynamics of a 
solid -propellant rocket; included are stress-strain laws for solid-propellant materials. 
Several problems are considered using a continuous representation of the solid- 
propellant rocket. 

Once a model for the solid-propellant rocket is obtained, the remainder of the 
vehicle can be modeled in the same manner as for a liquid-propellant vehicle. 

3 . 5 CLUSTERED BOOSTERS 

One method for obtaining the higher thrust required for large payloads is to attach 
rocket engines or motors to a central core; for liquid-propellant boosters a peripheral 
ring of propellant tanks is attached to a central tank and the engines are supported on 
truss members connecting the tanks; for solid-propellant boosters, the motors are at- 
tached to a central solid- or liquid -propellant booster. These clustered tank designs 
destroy axial symmetry and quite often planes of symmetry. Such configurations re- 
sult in a more complicated dynamic model where a number of cylindrical tanks are 
coupled by their elastic connections and must be allowed freedom in several directions 
for adequate description of vehicle modes (see Figure 24). 

For preliminary design it is sufficient to choose approximate planes of symmetry 
and analyze the vehicle for bending modes in pitch and yaw planes using branch beams 
connected to the central core by translational and rotational springs. Simplified tor- 
sional and longitudinal models will also suffice at this stage. These simple models 
can be used to identify possible problem areas (such as relative modal frequencies) 
and provide design criteria for the connections between tanks. 

A complete analysis (or carefully conducted test) should be undertaken to describe 
all the primary modes of the clustered vehicle. This analysis would provide transla- 
tion and rotation in two mutually pe rpend icular planes; torsion and longitudinal motion. 
The model of the tanks for translation and rotation in each of the two planes would be 
very similar to that discussed for the cylindrical booster. Provision must be made to 
account for the motion of the outer tanks in these two directions due to the torsional 
displacement of the central tank and the elastic connections. It is also possible that 
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longitudinal motion will couple with lateral 
and torsional displacement. As an example, 
consider a cluster arrangement where the 
connection at the bottom provides moment, 
shear, and axial restraints while the con- 
nection at the top provides only shear res- 
traint. Then it is possible to find a mode 
where the external tanks are bending, caus- 
ing moments and deflections at the connection 
to the center core which will result in longi- 
tudinal motion of the core. The significance 
of these types of modes can be ascertained 
only from the analysis (or test) and can vary 
greatly from vehicle to vehicle. 

The torsional properties in the model can 
be represented by the torsional stiffness and 
roll inertia of each tank. The tanks must 
then be connected by the elastic properties 
of the truss. The complete model for the 
clustered booster then consists of the axial 
load, shear, bending moment, and torque. 

The top connection transmits only shear. Be- 
cause of the nature of the connections, it can 
be seen that yaw bending and longitudinal 
coupling can occur. Pitch bending and tor- 
sion represent another possible coupling 
mechanism. Storey in Reference 14 develops 
the coupled flexibility matrices for these two 
conditions. This method encountered diffi- 
culty in that the number of stations required 
for adequate representation of the system 
with the required transformations exceeded 
computer capacity. 

The final Titan IIIC analysis presented 
in Reference 15 utilizes the mode synthesis 
approach. The longitudinal, torsional, and 
pitch and yaw bending modes are determined 
for each tank and are then coupled by the 
elasticity of the connecting elements. The 
influence coefficients for these trusses were 
obtained experimentally. A comparison of 
analytical and l/5-scale experimental re- 
sults is given in Reference 17. 
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3.6 LATERAL-TORSIONAL-LONGITUDINAL COUPLING 

The typical axisymmetric-cylindrical space vehicle is analyzed as if lateral, torsional, 
and longitudinal motion are not coupled. Actually, these vehicles are not completely 
symmetric and a possible coupling mechanism, however slight, can always be found. 
The importance of this coupling can vary greatly from vehicle to vehicle and even if it 
is known to exist from flight or experimental data, the coupling mechanism is difficult 
to identify. These coupling problems often occur when the modal frequencies of two 
modes, say, one lateral and one torsional, are very close together. Then a very small 
coupling mechanism, such as center-of-gravity offset from the supposed line of sym- 
metry, can result in coupled motion. 

A comparison of the frequencies of the modes in the three directions should be 
made to determine the existence of modes of nearly equal frequency. If such a condi- 
tion exists, it is necessary to examine the condition under which this may cause a sig- 
nificant problem. As an example, if excitation of a bending mode by an atmospheric 
disturbance occurred, could this cause excitation of a critical torsional mode at this 
same frequency ? 

Cylindrical vehicles with unsymmetric upper stages or payloads of large mass can 
cause coupling in the various directions in the low-frequency modes. The model and 
analysis then become complicated and approach that of the clustered boosters. Repre- 
sentation of this configuration requires detailed description in the unsymmetric stages 
and an analysis as described later for clustered boosters. Preliminary work would 
indicate the degree of sophistication to be used for adequate representation for loads 
analysis. 

The Saturn I vehicle consists of a center LO2 tank with eight peripheral tanks for 
alternating LC> 2 and RP-1. These tanks are connected at top and bottom by trusses 
providing axial, shear, and torsion restraint in both planes at the bottom plus moment 
restraint in the tangential planes. The top connection provides similar restraint ex- 
cept for the fuel tanks which do not transmit axial load. The trusses are not symmetric 
with respect to planes of symmetry of the tanks, but this effect is small so that planes 
of symmetry as defined by the tanks do not introduce large errors. 

Milner (Reference 18) establishes theoretically the uncoupling of pitch, yaw, and 
torsion modes for a symmetrical cluster booster and investigates the effect of minor 
asymmetry. Results of this study indicate that the effect of such coupling on natural 
frequencies is minor; mode shapes are not presented. 

Lianis (Reference 19) develops a matrix solution of the dynamics problem of a 
four-tank booster without center core. The flexibility matrix of the whole unit, with 
appropriate beam end fixity, is derived. This flexibility matrix together with a suit- 
able mass matrix is used to derive equations of free vibration in matrix form. The 
tanks are assumed to be similar, but the solution can be modified accordingly for the 
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case of nonsimilar tanks and for other tank configurations. The formulation is general 
so as to furnish any complex mode of vibration. Simple modes, however, can be ob- 
tained as particular cases of the general problem. 

3.7 DAMPING EFFECTS 

Dissipative (damping) forces exist in the vibrating structure as a result of material 
strain hysteresis and coulomb friction in structural joints. The nature of these damp- 
ing effects is obscure and does not lend itself to analysis other than an approximate 
empirical treatment, by which the gross effect of these scattered dissipative mecha- 
nisms is represented as equivalent viscous damping, added to each mode as appropri- 
ate. The damping is thus assumed to produce no coupling between modes. While this 
mechanization is not entirely realistic, it is justified by two observations: First, the 
actual damping is very low and is found by test to produce little coupling. Thus nearly 
pure normal modes of a system may be excited and the system observed to decay al- 
most harmonically. The indication given is that velocity-dependent coupling is very 
small. Second, if an attempt is made to show a velocity-dependent coupling, the co- 
efficient would have to be determined experimentally. Since the direct damping coef- 
ficient is itself difficult enough to measure it is clear that the accuracy of a study 
cannot be increased by the introduction of still more suspect data. The structural 
damping force is a function of the deflection of the generalized coordinate of the mode 
but in phase with the velocity of the generalized coordinate of that mode. To treat this 
damping as viscous damping requires that the mode oscillate in a quasi-harmonic man- 
ner. This damping force may then be expressed as a damping factor, where 
2C n co n | n is the internal damping force of the n^ mode per unit generalized mass. 

Fluid propellant damping forces result from the dissipative nature of a viscous 
fluid undergoing shear. Although there are some approximate methods for calculating 
damping forces, these forces are most commonly arrived at by testing the actual tank, 
in the case of small vehicles, and a model tank in the case of large vehicles. These 
forces may be represented as a propellant damping factor, C n » in the expression 
2C n to n 4 n which is the damping force per unit sloshing mass and | n is the lateral ve- 
locity of the n*"* 1 sloshing mass. 

Aerodynamic damping forces result from lateral velocity of the vehicle which 
causes, for any particular point on the vehicle, a small angle of attack. The aerody- 
namic force associated with this angle of attack opposes the lateral motion, thereby 
dissipating energy. The aerodynamic damping forces are easily calculated and are, of 
course, a function of 4 n . Aerodynamic damping on launch vehicles is often not impor- 
tant for dynamic load analyses; however, for some configurations, e.g. , hammerhead 
payloads and winged payloads, it may need to be considered. 
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4/METHODS for solution 


4.1 FORMATION OF THE EQUATIONS FOR SOLUTION 

The longitudinal model provides a mathematical representation of the real physical 
system. As outlined in Section 3, dynamic characteristics of the model are deter- 
mined from the mass, stiffness, and dynamic matrices. 

4.1.1 STIFFNESS MATRIX. Formation of the dynamic matrix for computation of 
modal properties requires formation of the mass and flexibility matrices. The flexi- 
bility matrix is not ordinarily formed directly but is usually obtained by inversion of 
the stiffness matrix. (See the Reference 8 matrix inversion methods.) This proce- 
dure is followed since the stiffness matrix can generally be formed in a simpler and 
more direct manner than the flexibility matrix (see Reference 6). 


Where spring-mass models are employ- 
ed, the relationship between spring constants 
and the stiffness matrix is fairly simple. 

The stiffness matrix element K^j is the force 
at i due to a unit displacement at j only with 
all other displacements equal to zero. Off- 
diagonal elements K^j of the stiffness matrix 
are simply the negative value of the stiffness 
of the spring connecting node i to the node j 
of the model. The diagonal element of 
the stiffness matrix represents the sum of 
the stiffnesses of all springs connected to 
node i. As an example, consider the spring- 
mass model of Figure 25. 



The stiffness matrix for the model re- Figure 25> Spring _ Mass Mode i 

lates forces to displacements as follows. 
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F 3 


0 

- K 23 

K 23 + K 34 

~ K 34 


X 3 

s 

1 

0 

- K 24 

- K 34 

K 24 + K 34_ 


*4. 


Note that the coefficients in column one are forces due to a unit displacement at node 
1 only with displacements at 2, 3, and 4 of zero. The force required to produce the 
unit displacement at node 1 is while -K^ is the reaction at 2 corresponding to the 
unit displacement at 1. 
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If the following stiffness matrix is for absolute displacements of a free structure, 


[K] 


K 11 

K 12 

* 

CO 

K 21 

K 22 

K 23 

K 31 

K 32 

K 33 


(118) 


and the associated mass matrix is a diagonal matrix, the corresponding model and 
springs are 


K 

-K 

12 

12 

K 

-K 

13 

13 

K 

-K 

23 

23 
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where the model is given in Figure 26, 

The methods for deriving stiffness matrices or 
spring-mass models for propellant tanks and other por- 
tions of the structure were discussed in Section 3. Assem- 
bly of the vehicle stiffness matrix is accomplished by su- 
perposing or adding the stiffness matrices of the individual 
elements . 


'12 


‘13 


'23 


Figure 26. Model Corres- 
ponding to 
Stiffness Matrix 


One of the problems of obtaining a solution to the 
equations of motion for a free-free structure with no 
external forces applied is that the stiffness matrix of 
such a structure is singular. The structure could move 
as a rigid body upon application of any external force. 

In this instance inverting the stiffness matrix to obtain a flexibility matrix may not be 
done directly. Several methods are available to overcome this difficulty. They are 
essentially the use of an additional equation, which states that for free vibration of a 
free-free structure the summation of the inertial forces is zero: 


m 


- co 


£ m i x i 


= 0 


(120) 


i=l 


From this equation we may choose an arbitrary mass, say m^, and let 
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(121) 



Then in matrix notation, 

X 1 = - m i’ 1<m 2 m 3 ••••'• m m > {X} 

= <E> fx'} 

fx} = 

Then 



Then from Equation 12 , page 7 , with the stiffness given by 



Expanding the second equation , 

-to 2 [m'] {x} +|{y} <E> +[/3]]tx'} = 0 

Therefore 

[K'] = {y} <E> + m 



(122) 

(123) 

(124) 

(125) 

(126) 

(127) 

(128) 
(129) 
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and 


—2 fa*) = t K '] 1 [ m 'l fV] (130) 

co 

The above equation may be used to solve for the natural modes and frequencies of 
the free structure. Once x' has been determined, may be determined from 
Equation 123. 

4. 1.2 FLEXIBILITY MATRIX. The flexibility matrix may be developed directly for 
simple statically determinant structures by determining the displacements due to unit 
loads. However, obtaining the stiffness matrix and inverting it is easier when very 
many redundancies exist in the structure. For direct development of the flexibility 
matrix in terms of the flexibility matrices of individual structural elements for com- 
plex structures, see Reference 20. 

4.1.3 TRANSFORMED MASS MATRIX. An approach that is particularly advantageous 
tor close-coupled systems is that of transforming the coordinate system from the abso- 
lute to the relative sense. In prior discussions the displacements of the system co- 
ordinates have been referenced to a fixed point or neutral position. These same dis- 
placements may also be expressed relatively; referenced to an adjacent coordinate. 

The relationship between displacements in absolute terms, {x} , and the displacements 
in relative terms, [x] , is readily described by a simple transformation matrix: 

{x} = [T] { x } ( 131 ) 

Thus the kinetic energy of the system can be expressed in terms of relative co- 
ordinates by 


2KE = {x}* [T]' [M] [T] fx] (132) 

The deflections of the connecting springs are expressed in terms of relative co- 
ordinates, which in its most general form requires another transformation matrix. 

{A} = [TR] [x] (133) 

If the number of springs is equal to or less than the number of inertias, the de- 
flection of each spring will be defined by a different relative displacement; consequently 
the transformation matrix, [TR], can be written as a diagonal matrix of unit elements. 
For such a case, [TR] may be neglected without affecting the solution. 


The potential energy of the system is given by 
2 PE = {5}'[TR]' [K] [TR] {£} 
where K is a diagonal matrix of the spring rates. 


(134) 
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Substituting the kinetic energy and potential energy terms into LaGrange's equation 


d [ a(KE) l + apE + aw 
dt [ 34. dq. 


(135) 


the equations of motion become 

[T]' [M] [ T] { 3c } + [TR] ' [K] [TR] {*} = 0 (136) 

and the dynamic matrix is established by 

[([TR]' [K] [TR]) _1 [T]' [M] [T] ] {0} = A {0} ( 137) 

The advantages of this method are greatest when the system is free or cantilever- 
ed and multiple load paths are absent or constitute only a minor portion of the system. 
Under these conditions, the transform matrix [TR] contains little or no off-diagonal 
terms; the matrix [TR] ' [K] [TR] is no more than slightly coupled and may be in- 
verted with a minimum of effort. 

When the number of springs is equal to or less than the number of inertias, the 
transform matrix is a diagonal matrix of unit elements; [TR] [K] [TR] reduces to 
[K] which, being a diagonal matrix, can be immediately inverted by taking the reci- 
procals of the individual elements of the diagonal. The dynamic matrix can be deter- 
mined with a minimum of effort and operated on to obtain modal data. 

4.2 SOLUTION FOR CHARACTERISTICS 

After the stiffness and mass matrices for the model are formed, the equations of 
motion may be written with external forces equal to zero. This set of linear differ- 
ential equations may be solved for modal properties by several techniques. If the 
equations for harmonic motion , 

-u 2 [M] [x] +[K] [x] = 0 (138) 

are written in the form 

|[K] - oj 2 [M]j [x] = 0 ( 139 

then for 

[x] ft 0 
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a solution exists only if the determinant 


[K] - CO 2 [M] 


0 


(140) 


Expansion of the above determinant leads to a polynomial. If [K] and [M] are n^ 1 
order matrices , the polynomial (characteristic equation) will be n^ 1 order and have 
n roots. The roots of the polynomial yield the natural frequencies u). Particular 
natural frequencies may be inserted in Equation 139 to determine the corresponding 
natural mode shape. This procedure results in determination of all natural frequencies 
and mode shapes. The method is very useful for small matrices but becomes difficult 
to use for large matrices. Since only lower frequency modes are generally of impor- 
tance, other methods are usually used which do not require direct solution for the roots 
of the characteristic equation. Among these methods are the matrix iteration, Holtzer 
(Myklestad), and energy (Rayleigh-Ritz) methods. Matrix iteration is discussed in 
References 8 and 21; the Holzer (Myklestad) method and the Rayleigh-Ritz method are 
discussed in References 21 and 22. 


4.3 MODAL SYNTHESIS 


If a structural model requires a very large number of degrees of freedom, compu- 
tation of modal properties of the entire structure directly may pose a difficult pro- 
blem. The structure could be divided into several parts and the modal properties 
could be computed for several modes of each part. The dynamic properties of the 
entire structure could then be determined from the properties of the parts using the 
method of modal synthesis. This method is discussed in detail in References 13, 21, 
and 23. 


56 



5/REFERENCES 


1. L. L. Fontenot and G. Lianis, "The Free Vibrations of Thin Elastic Pressurized 
Cylindrical Shells with a Free and Incompressible Liquid Having a Free Surface, " 
International Symposium on Space Technology and Science , Paper No. 6-6-8, 1963. 

2. J. D. Wood, Survey on Missile Structural Dynamics , Report 7 102-004 1-NU -000, 
Contract No. AF 04(647)-619, 1 June 1961. 

3. R. G. Rose, J. A. Staley, and A. K. Simson, A Study of System-Coupled Longi- 
tudinal Instabilities in Liquid Rockets, Report AFRPL-TR -65-163, Contract AF 
04(611)-9956, September 1965. 

4. L. D. Pinson, Longitudinal Spring Constants for Liquid-Propellant Tanks with 
Ellipsoidal Ends , NASA Technical Note D-2220, November 1964. 

5. J. S. Archer and C. P. Rubin, Improved Analytic Longitudinal Response Analysis 
for Axisymmetric Launch Vehicles, NASA CR-345, Contract NAS1-4351, 

December 1965. 

6. M. J. Turner, R. W. Clough, H. C. Martin, and L. J. Topp, "Stiffness and 
Deflection Analysis of Complex Structures," Journal of the Aeronautical Sciences , 
Vol. 23, No. 9, September 1956. 

7. J. H. Baltrukonis, A Survey of Structural Dynamics of Solid Propellant Rocket 
Motors , NASA Report NASA CR-658, December 1966. 

8. V. N. Faddeeva, Computational Methods in Linear Algebra , Dover Publications, 
New York, 1959. 

9. S. Timoshenko, Theory of Plates and Shells, McGraw-Hill Book Company, Inc., 
First Edition, 1940. 

10. R. G. Rose, Dynamics of the Atlas 5-CPS Lo ngitudinal Oscillation Following 
Launch as Related to the Tank Pressure Regulation System , Convair Division 
Report GDA63-0712, Volumes I and n, Contracts AF 04(694)-196 and AF 04(695)- 
240. 

11. L. D. Pinson, H. W. Leonard, and J. P. Ramey, "Analysis of Longitudinal 
Dynamics of Launch Vehicles with Application to a l/10-Scale Saturn V Model," 
AIAA/ASME 8th Structures, Structural Dynamics, and Materials Conference, 

Palm Springs, March 1967. 


57 





12. R. H. Schuett, B. A. Appleby, J. D. Martin, Dynamic Loads Analysis of Space 
Vehicle Systems, Launch and Exit Phase, Convair division of General Dynamics 
Report GDC-DDE66-012, June 1966. 

13. W. C. Hurty, Dynamic Analysis of Structural Systems by Component Mode 
Synthesis , TR 32-530, Jet Propulsion Laboratory, Pasadena, California, 
January 1964. 

14. S. I. Gravitz, "An Analytical Procedure for Orthogonalization of Experimentally 
Measured Modes, " Journal of the Aerospace Sciences , November 1958. 

15. R. E. Storey, "Dynamic Analysis of Clustered Boosters with Application to 
Titan III," AIAA 1963 Summer Meeting, Paper No. 63-208, June 1963. 

16. Bodley, Ikard, and Schulta, Vibration Analysis Report, Program 624A, Con- 
figuration C, Flight Plan VIP , Air Force Report No. SSD-CR-65-1, January 
196b. 

17. I. J. Jaszlics, Program 624A 20%-Scale Dynamic Model Ground Vibration Sur- 
vey , Final Test Report, Vol. Ill, SSD-CR -65-47, Martin Company, Denver, 
Colorado, May 1965. 

18. J. L. Milner, Three Dimensional Multiple Beam Analysis of a Saturn 1 Vehicle, 
NASA TMX-53098, July 1964. 

19. G. Lianis, Matrix Analysis of Vibrations of Clustered Boosters , Convair 
Division Report AE6 1-0858, September 1961. 

20. E. C. Pestal and F. A. Leckie, Matrix Methods in Elastomechanics , McGraw- 
Hill, 1963. 

21. A Monograph on Lateral Vibration Modes , Convair Report GDC -DDF65-001, 
Contract NAS8-11486, 22 February 1965. 

22. L. S. Jacobsen and R. S. Ayre, Engineering Vibrations , McGraw-Hill, 1958. 

23. A Monograph on Torsional Vibration Modes , Convair Report GDC-DDF65-003, 
Contract NAS8-11486, 3 May 1965. 


58 


NASA-Langley, 1967 32 CR-9J6 








